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Abstract 

We use numerical techniques to study the formation of singularities 
in Ricci flow. Comparing the Ricci flows corresponding to a one pa- 
rameter family of initial geometries on 5 3 with varying amounts of S 2 
neck pinching, we find critical behavior at the threshold of singularity 
formation 



1 Introduction 

Given a manifold M n and a Riemannian metric g specified on M n , the Ricci 
flow determines a one parameter family of metrics g(t) via the geometric 
evolution equation 

d t g ab (t) = -2R ab [g{t)l (1) 

with initial condition 

9(0) =9- (2) 

*Email: garfinkl@oakland.edu 
'Email: jim@newton.uoregon.edu 



1 



(Here is the Ricci tensor of the metric g a b.) Ricci flow has been a very 
effective tool for studying the sorts of special geometries which a manifold 
admits. In most of these applications to date [S| [7], the flow is shown 
to converge to the special geometries (e.g., to a constant negative curvature 
hyperbolic metric for the Ricci flow on any two dimensional manifold of genus 
greater than one). Future applications, however, are likely to require that 
one understand Ricci flows which develop singularities. 

In this work, we begin a program of numerical study of Ricci flow singular- 
ities. We focus here on the following question: Say we have a one parameter 
family of metrics g\ specified on the manifold S 3 and say we know that the 
Ricci flows g\(t) starting at g\ for very large A converge (with the volume 
suitably normalized) to the round sphere metric, while the Ricci flows g\(t) 
for small values of the parameter A become singular, in the sense that (re- 
gardless of volume norm) the curvature of g\ (t) grows without bound as t 
increases. What happens to the Ricci flows g\{t) with intermediate values 
of A? Is there, in particular, a certain threshold value \ cr it for which the 
Ricci flow has interesting intermediate behavior: neither convergence, nor 
formation of a standard singularity? 

For solutions of Einstein's equations representing gravitational collapse, 
this sort of question has been studied extensively pQ 0j. Very distinct thresh- 
old behavior has been found, with a remarkable degree of universality. That 
is, in examining a number of one parameter families of initial data for gravita- 
tional collapse, one finds qualitatively the same sort of discretely self-similar 
solution occurring for threshold initial data for all of these families. 

For Ricci flow, our work here is the first search for critical behavior. We 
have examined a particular one parameter family of spherically symmetric 
"corseted sphere" geometries on S 3 , with the parameter describing the degree 
of corseting at the equator, and therefore parametrizing the magnitude of the 
S 2 neck pinch curvature at the equator. We do indeed find numerically that 
for geometries with a small amount of corseting, the Ricci flow converges to 
the round sphere metric, while for geometries with a large amount of corset- 
ing, an S 2 neck pinch singularity occurs. Moreover, we find that there is a 
critical value of the parameter, dividing the two regimes. Finally, our studies 
show that the Ricci flow for the geometry marked by the critical parame- 
ter value neither converges to the round sphere geometry nor forms an S 2 
neck pinch singularity. Instead, this flow approaches a "javelin" geometry, 
marked by curvature singularities at the poles, with roughly uniform curva- 
ture between the poles. This javelin geometry corresponds to the "type 3" 
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singularity described by Hamilton [8] and discussed by Chow. [2] 

We describe in detail in Section 2 the corseted sphere geometries that we 
study here, and write out the detailed form of the Ricci flow equations for 
these geometries. Also in Section 2 we describe our numerical methods. We 
present our results in Section 3, noting the behavior of the Ricci flows for 
subcritical, supercritical, and critical initial geometries. Concluding remarks 
appear in Section 4. 

2 Corseted Sphere Geometries 
and Their Flow Equations 

The corseted sphere geometries and their flows are all represented by spher- 
ically symmetric metrics on S 3 of the form 

g = e 2X [e~ 2W difj 2 + e 2W sm 2 if)[d6 2 + sinW]) (3) 

Here (if), 8, <f>) are standard angular coordinates on the three sphere; spherical 
symmetry holds so long as we assume that the metric functions X and W 
are functions only of if). 

Smoothness of the metric at the poles, where if) takes the values and 
7r, requires that X and W be even functions of sinif) in a neighborhood of 
the poles. Therefore d^X and d^W must vanish at the poles. However, 
smoothness of the metric also requires that W vanish at the poles. To avoid 
the numerically awkward imposition of two conditions on W at a single point, 
it is convenient to replace the variable W by S = W/siv^ip. Smoothness of 
the metric at the poles can then be enforced by the requirement that d^X 
and d^S vanish at if) = and if) = it. 

To obtain the corseted sphere geometries, we set W = X, and choose X so 
that 4e 4X sin 2 ^ = sm 2 2ip for cos 2 '?/' > 1/2 and 4e 4X sin 2 ?/> = sin 2 2?/> + 4Acos 2 2'0 
for cos 2 ip < 1/2. Here A is a constant, which parametrizes the degree of 
corseting for these geometries. For A = 0, the geometry represents two round 
three spheres joined at the poles. This is a singular geometry. For A positive, 
the cusp smooths out and the geometry is non singular; however, for small 
values of A one still expects the curvature at the tightly pinched equator to be 
very large. In Figure 1, we graph the area of the if) = constant cross-sections 
as a function of if> for a few representative values of A. In Figure 2, we 
graph the Ricci curvature eigenvalue in the direction along the S' 2 symmetry 
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Figure 1: area of the S 2 s as a function of ip for A = 0.1 (solid line), A = 0.15 
(dashed line) and A = 0.2 (dotted line). 
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Figure 2: R S 2 as a function of ip for A = 0.1 (solid line), A = 0.15 (dashed 
line) and A = 0.2 (dotted line). 

as a function of ip. Note that as a consequence of the assumption in these 
geometries that W = X, we verify that the coordinate ip gives the value of 
the geodesic distance in the radial direction. 

We are interested in the Ricci flow of the corseted sphere metrics. How- 
ever, since the Ricci flow equation is only weakly parabolic, and since 
numerical evolutions appear to be more stable for strongly parabolic sys- 
tems, we instead work with the DeTurck flow. [3 3 The strongly parabolic 
PDE system generating the DeTurck flow is 

d t g ab = -2R ab + 2D (a V b) (4) 

where D a is the derivative operator associated with the metric g ab , and where 
the vector field V a is given by 

V a = g bc (f% c - A a bc ) (5) 

with being the connection of the metric g ab and with A£ c being any fixed 
connection. For a given initial geometry, one can find the corresponding 
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Ricci flow g(t) by first finding the corresponding DeTurck flow g a b, and then 
pulling back via the time dependent diffeomorphism generated by the vector 
field V defined above. 

Note that neither the Ricci flow ([TJ nor the DeTurck flow (@J) preserves 
volume. One can normalize the volume for either of the flows by adding 
the term y g a b, with f the spatial average of the scalar curvature, to the 
flow equation. Alternatively, one can control the volume along either flow by 
periodic uniform blowups. 

We now calculate the volume normalized DeTurck flow for spherically 
symmetric metrics of the form (jSJ). Choosing the reference connection to be 
that of the round sphere, and using primes to denote the spatial derivative 
d^p, we obtain the evolution equations 



d f X = e 2 ^"^ 



X" + 2cot^X' - 2 + -([Xf + [W] 2 ) + 3X'W 



+ (1 - e~ m ) (— Vr + l + 2cotVW 
V 2sm ip , 



r 

+ 3. (6) 



W" + 2 cot xj}W' - ^{[Xf + [Wf) - 3X' W 



(7) 



+(1 - e- AW ) ( 1 - 2 cot MX' + 2W'\ 

Here we note that the average scalar curvature f is given by 

f = — f r di[>e x+3W (e~ m - 1 - 4 sin ^ cos i[>W' + sin 2 ^[3 + (X' + W') 2 ]) , 
N Jo ' 

(8) 

where the normalization constant N is given by 

N= f dip e 3X+w sin 2 ?/*. (9) 
Jo 

We also note that the only non vanishing component of the vector field V 
defined in © is 

V i , = -l(3W + X' + 2cotiP[l-e- iW }) (10) 

As discussed earlier, it is useful to work with the quantity S = W/sm 2 ip, 
rather than W. Since this definition implies that W = sm 2 ipS'+2 sin ip cos ipS, 
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we readily obtain 

dt S = e 2 ^"*) 



1 - e~ 4W 



S" + 6cot^' -85 ^- fl -AW -e~ m 

2sm ip v 

1 - 2 [cot ij)X' + 2 sin^ cos ^5' + 4cos 2 ^>S] 



sin 2 -?/' 

/ [X'/ sin V] 2 + [sin + 2 cosVS'] 2 + 6[X'/ sin -0] [sin + 2 cos 



.(11) 



Equations (jOJ and are the ones that we evolve numerically, noting that 
wherever we encounter W or d^W these are to be expressed in terms of S 
and d^S. 

3 Numerical Results 

To study these evolution equations numerically, we proceed as follows: We 
divide the spatial coordinate range (0, ir) of ip into N — 2 pieces, so that we 
have Aip = it/ (N — 2). We choose iV grid points, including a pair which run 
outside the coordinate range. So the first spatial grid point is at tp = —Aip/2, 
while the last is at ip = n + Aip/2. Then a function of the form F(ip, to) for 
fixed time t is replaced by a set of N numbers = F((i — |) Atp, t ) where 
1 < i < N. Spatial derivatives are replaced by centered finite differences in 
the following way: 

d^F -> ——2 (13) 

For the time dependence of these functions, we choose a fixed time step At, 
and replace F((i — ¥)Aip, nAt) by the numbers F™. 

Now, for an evolution equation of the form dtF(tp,t) = G{ip,t) we nu- 
merically evolve using the approximation 

F n+1 =F n + AfG n (14) 

This evolution is implemented for all values of i except 1 and N. Note that 
these two "ghost zones" are not part of the manifold since ip is not in the 
range < ip < n. At the ghost zones we use smoothness of the metric which 
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Figure 3: X for subcritical Ricci flow 



implies that X' and S' vanish at the poles. We implement this condition as 
X 1 = X 2 and X N = X N _i (and correspondingly for S). 

Runs were done (on Unix and Linux workstations), starting from a wide 
range of initial corseted sphere geometries (parametrized by A). As noted 
above, for small A we expect that the geometry is sufficiently close to singular 
that the distorted S 3 will pinch off into a singularity, while for sufficiently 
large A the flow may overcome the distortion to evolve the data to a single 
round S 3 . 

These expectations are confirmed by our numerical simulations. Figures 
3 and 4 show the results of a "subcritical" run, i. e. one which does not result 
in a singularity. For this run we have chosen A = 0.2. Note that at late times 
we find that X approaches a constant and that S — > 0. These are the values 
for a round S* 3 . 

Rather than focusing on the evolution of the metric components, it is 
more instructive to examine the behavior of curvature. As a consequence 
of the spherical symmetry of the geometries, the Ricci tensor R a i, has two 
independent eigenvalues. We call these eigenvalues R S 2 and R± where R S 2 
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corresponds to the eigenspace in the symmetry S 2 directions and R± corre- 
sponds to the eigenspace orthogonal to the symmetry S 2 directions. Some 
straightforward calculations show 



R ± = -2e 2(w - x ^ [-1 + X" + W" + (X' + 3W') cot V + 2(X' + W')W'\ (15) 

R s2 = - e m~X) _ 2 + + X n + w „ + , zx , + 5W n CQt ^ 

sin ip 

+(X' + W')(X' + 3W') (16) 

(where we note that W, W and W" can be expressed in terms of S, S' and 
S"). One can express the invariants of the Ricci tensor (and since we are in 
3 dimensions, the invariants of the Riemann tensor as well) in terms of Rsz 
and R±. We have 

R = 2R S2 + R ± (17) 
R ab R a b = 2R 2 S 2 + R\ (18) 
R abcd Rabcd = 4R ab R ab -R 2 = 2R 2 ± + (R ± - 2R S2 ) 2 (19) 

Figures 5 and 6 show the behavior of R s i and R± for the same subcritical run. 
Note that both eigenvalues asymptotically approach the same constant value 
at late times, which confirms the contention that the flow of a subcritical 
geometry converges to the round sphere geometry. 

We next consider the DeTurck flow of "supercritical" initial geometries, 
i.e. those for which the evolution is singular. Figures 7 and 8 plot the 
behavior of Rs^ and R± for an initial geometry with A = 0.11. Characteris- 
tically, we find that in the neighborhood of the equator, i?^2 grows without 
bound as t increases, while R± appears to stay bounded. As well, both are 
bounded away from the equator. This signals the formation of an S 2 neck 
pinch singularity at the equator. 

The behavior just described for the DeTurck flow-and consequently for 
the Ricci flow-for A = 0.11 corseted sphere data is also found for any initial 
geometry with A < 0.11. Similarly, for A > 0.2, the flow has the subcritical 
behavior illustrated in Figures 5 and 6. This has led us to seek critical 
behavior at a threshold value. Using a binary search, we have located the 
threshold value at approximately A = 0.1639 

The DeTurck flow starting at a corseted sphere geometry with this value 
of A behaves differently than the flow for both subcritical and supercritical 



10 



Figure 5: R S 2 for subcritical Ricci flow 
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Figure 6: R± for subcritical Ricci flow 
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Figure 7: _R 5 2 for supercritical Ricci flow 
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Figure 8: R± for supercritical Ricci flow 
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Figure 9: R S 2 for critical Ricci flow 
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Figure 10: R± for critical Ricci flow 
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geometries. As seen in Figures 9 and 10, R±_ gets small everywhere except 
at the poles, while Rg^ slowly grows every where except at the poles. At the 
poles, both curvatures get very large. In a sense, the geometry approaches 
that of a three dimensional javelin, with S 2 cross-sections. 

4 Conclusions 

The numerical work we have described here clearly shows rather special be- 
havior for the Ricci flow at the threshold parameter value for a one parameter 
family of corseted spheres. Is this behavior in any sense universal? This has 
not yet been determined. We plan to examine other families of initial ge- 
ometries to see if the threshold behavior persists. One interesting set of 
geometries we plan to examine are those with, initially, more than one neck 
pinch.. Do the neck pinches coalesce? Do we get javelin geometries for 
threshold initial data? We also plan to consider families of initial data which 
are not spherically symmetric. 
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